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Abstract 

We study the Type OA string theory in the (2, 4/c) superconformal minimal model back- 
grounds, focusing on the fully non-perturbative string equations which define the partition 
function of the model. The equations admit a parameter, T, which in the spacetime in- 
terpretation controls the number of background D-branes, or R-R flux units, depending 
upon which weak coupling regime is taken. We study the properties of the string equa- 
tions (often focusing on the (2,4) model in particular) and their physical solutions. The 
solutions are the potential for an associated Schrodinger problem whose wavefunction is 
that of an extended D-brane probe. We perform a numerical study of the spectrum of 
this system for varying T and establish that when T is a positive integer the equations' 
solutions have special properties consistent with the spacetime interpretation. We also 
show that a natural solution-generating transformation (that changes F by an integer) is 
the Backlund transformation of the KdV hierarchy specialized to (scale invariant) solitons 
at zero velocity. Our results suggest that the localized D~branes of the minimal string 
theories are directly related to the solitons of the KdV hierarchy. Further, we observe an 
interesting transition when F = — 1. 



^Also, Visiting Professor at the Centre for Particle Theory''. 



1 Background 



The (2,4/c) series of the minimal type OA string theory, in the presence of background R-R 
sources, has a non-perturbative definition via the following "string equation" pQ : 

M7^2-i7^7^" + i(7^')2 = ^.2^2 . (i) 

Here u{z) is a real function of the real variable z; a prime denotes ud/dz; and F and u are 
constants. The quantity 71 is defined by: 



k=0 ^ ^ 



(2) 



where the Rk {k = 0, . . .) are polynomials in u{z) and its 2;-derivatives. They are related by a 
recursion relation: 

-Rfe+i = ~ - -ju'Rk , (3) 

and are fixed by the constant Rq, and the requirement that the rest vanish for vanishing u. 
The first few are: 

Ro = \; Ri = -\u- R2 = 4(3w' - u) . (4) 
2 4 lb 

The kth. model is chosen by setting all the other ts to zero except to = z, and tk, the latter 
being fixed to a numerical value such that TZ = T>k — z. The Vk are normalised such that the 
coefficient of is unity, e.g.: 

Pi = M , Do = u u , V-x = u — uu iu) H u . (5) 

' 3 2^ ^ 10 ^ ' 

For the fcth model, equation has asymptotics: 

/ \ 1 '^F 

= zk ^ r + "' ^ — > +cxd , 

z/2(4F2 - 1) 

uiz) = h ■ ■ ■ for z — > —00 . (6) 

Az^ 

The function u{z) defines the partition function Z = exp(— F) of the string theory via: 



u{z) = V 



(7) 



where is the coefficient of the lowest dimension operator in the world-sheet theory. Integrating 
twice, the asymptotic expansions in equations © furnish the partition function perturbatively 



2 



as an expansion in the dimensionless string coupling 

V 



9s 



i+— 



The role of F is now clear: at large positive z we find that F controls the number of background 
D-branes since the perturbative series contains both open and closed string worldsheets with 
a factor of F for each boundary; while the large negative z series gives only closed string 
worldsheets, with F^ appearing when there is an insertion of pure R-R fiuxmi2]. 

From the point of view of the kth theory, the other t^s represent coupling to closed string 
operators Ok- It is well known that the insertion of each operator can be expressed in terms of 
the KdV fiowsHHH: 

du , 

fife = • (8' 

The operator Oq couples to to? which is in fact — 4z, the cosmological constant (in the unitary 
model). So is often referred to as the puncture operator, which yields the area of a surface 
by fixing a point which is then integrated over in the path integral. So u{z) is the two-point 
function of the puncture operator. 

For F = the string equation (Q) was discovered by defining a family of string theories using 
double scaled models of a complex matrix M: 



Z = j dMexp!^-^TT{V{MM^)) 



(10) 



The resulting physics captured by the string equation was the first complete non-perturbative 
definition of a string theory|SJ|Hl[ZlElin]; sharing the large z perturbation theory of the original 
(bosonic) string theories obtained by double scaling of Hermitian matrix modelsjTIH ITU IT^ . 
but not suffering from their non-perturbative shortcomings. 

Non-zero F can also be studied using a matrix model definition in a variety of ways. One way is 
to add a logarithmic term|13| IT^ to an appropriate matrix model potential. Expanding 
the logarithm clearly adds holes of all sizes to the string worldsheets defined by studying the 
duals of the Feynman diagrams of the model. Another (equivalent) method is to define M as 
an (A^ + F) X rectangular matrix^ll2]. The double scaling limit (in which A^ is taken large) 
then yields equation ((T)). 

It is clear from these two methods that, in the interpretation of the matrix model as the world- 
volume theory of A^ D-branes (which defines the closed string theory holographicallv pTj after 
taking the double scaling limit), the modification corresponds to adding F "quark fiavours" to 
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the world-volume model. This corresponds to adding sectors of open strings stretching between 
the D-branes and F extra D-branes. 

So from this point of view it appears that F is a positive integer. However, this is not at all 
clear from the string equation itself. In fact, it is known perturbatively that the solutions of 
the equation have special properties for various fractional values of F ^. These might well turn 
out to be unphysical values, but this is not a priori clear. One of the purposes of this paper 
is to show that F being a positive integer is in fact selected out as special from the point of 
view of the string equation. We will do this by studying the equation non-perturbatively using 
numerical methods, and by also observing the properties of a natural transformation between 
solutions of the string equation which change F by ±1. This transformation was discovered in 
refpP, and we derive an explicit form for it here, and employ this form in our investigations. 
We further observe that the transformation is in fact a specialisation (to the scale invariant 
soliton sector) of the Backlund transformations of the KdV hierarchy! 

These Backlund transformations are known to change the soliton number of solutions of the 
KdV equation (and its higher order analogues) |18j. and so our observation is particularly in- 
triguing, since it implies that the (localized ^^1) D-branes of the minimal type OA string theory 
are intimately connected with the solitons of the KdV hierarchy. 

Note: While this paper was being written, a paper appeared|2ni which contains results that 
have some overlap with ours. 

2 A Spectral Problem 

A solution u{z) of the string equation serves as the potential for a one-dimensional Hamiltonian 
Ti which arises naturally in the double-scaled matrix modeipT]: 



In fact, the (first derivative of) the string equation can be obtained by eliminating the wave- 

^Note in particular the case F — ±1/2 in the large negative z regime of equation jnj. There are several 
other interesting cases of this type, including a family of exact double pole solutions, for half-integer V. See 





where v^z/ plays the role of h. 



refs.pCni. 
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function ip{z) from the equations: 




■s 



(12) 



Here, P is a differential operator representing scale transformations. It is constructed from (the 
differential operator part of) fractional powers [221 123 of Q, as follows 0: 



which yields an equation for n(z) which is the first derivative of equation 

The spectrum of this model has interesting and important information, and is a rather direct 
way of getting access to the non-perturbative physics of the modeP. The problem is naturally 
defined in terms of the eigenvalues. A, of the combination MM^. These are naturally positive. 
So in the double scaling limit (which focuses on the infinitesimal region near the tail of the 
eigenvalue distribution of the model) the scaled eigenvalues, which we denote as A^, are also 
distributed on the positive real line. 

We verified this in the full non-perturbative regime by studying the spectrum of 7^ numerically. 
To do this we solved the differential equation for a given value of /c as a boundary value problem, 
with the perturbative boundary conditions given in equations (jH)). We started out by using the 
equation-solving routine dsolve in the package Maple 9, although later we used other methods 
(see below for the case of negative F). We used a discrete z lattice of up to 8192 points for 
positive r. 

First, note the form of the function u{z) for some strictly positive values of F (see figure QJ. 
Since the potential tends to zero at ^ — — oo, and rises monotonically in the z — » +oo limit, it 
is clear that the spectrum has the chance to be both continuous and bounded from below by 
zero, as we expect from perturbation theory. A striking feature is the fact that the potential 
develops a potential well in the interior, as F ^ 0. 

It is therefore interesting and important to non-perturbatively establish these properties (conti- 
nuity, boundedness) of the spectrum for various F. First, we checked the spectrum of the case 

^The wavefunction i/j is a natural extended D-branej21I probe of the model, a fact which has been studied 
extensively recently for the bosonic string in ref. 25 . 




k=l 



(13) 



Eliminating ip, one has the operator equationj^j: 



[P,Q] = Q , 



(14) 
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Figure 1: Numerical solutions to equation for u{z): (a) The case fc = 1, and (b) the case k = 2, 
for a range of values of the parameter T. The very bottom curve in each case has F = 0. The solutions 
all asymptote to u = z^^^ for large positive z and u = for large negative z. 

r = (the case with no background branes), and verified that there are no bound states, despite 
the appearance of the well^^. We set up the problem numerically by discretizing the spatial 
coordinate typically using the interval —10 < z < +10, with u = 1. We used a range of 
lattice spacings in studying our problem to ensure stability. This was equivalent to a number of 
lattice points ranging from as little as 500 up to 10000 (the latter used when more accuracy was 
needed, as we shall see later on.) Discretizing this problem turns it into a problem of diagonal- 
izing a tri-diagonal matrix for which we coded a fairly efficient well-known method (the TQLI 
method in C++ to do the work. This allowed us to establish non-perturbatively with 
some confidence that there are no bound states in the well for positive F, and so the spectrum 
is continuous and bounded from below by zero, in accord with perturbative expectations. 

The next step is the case of negative F, which is clearly allowed by the equation and so should 
be studied. Our first observation is that solving the differential equation numerically became 
much more prone to error, and we had to use much more accurate methods. We used the NAG 

■^That this is the case exactly fits with a suggestive (but rough) estimate one can make by comparing the 
depth of the well with the value of the lowest energy eigenstate of the equivalent harmonic oscillator, determined 
by reading off the value of the second derivative at the bottom of the well (computed numerically). The well 
falls short in depth by ~ 10%. 

^This complements earlier studies of the properties of the F = solutions carried out extensively in refs-O 

El 13 El EI 
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library of routines for solving boundary value problems with C++ to make further progress, 
as they gave much greater control over the problem. We were able to construct directly in this 
way a solution for u{z) for values as low as F = —0.970 before running into numerical problems. 
See figure 121 




Figure 2: Numerical solutions to equation for u{z) in the case k = 1, for a range of negative values 
of the parameter F: 0.000, -0.500, -0.950, -0.965, -0.968, -0.970. 

Once again we should check that the perturbative expectation that the spectrum is positive 
persists non-perturbatively. Several solutions u{z), for various F down to —0.970 were used as 
potentials for Ti. and found (upon placing in our tridiagonal matrix diagonalization routine) to 
have no bound states. 

With these methods, the growing instability is suggestive, but not compellingly so, that some- 
thing interesting might be happening at F = — 1. In order to establish this more satisfyingly we 
employed some exact methods (in the form of a solution-generating transformation), which we 
will describe next. This will allow us to generate new u{z) from the ones we already can extract 
numerically, and allow us to approach F = — 1 in a controlled manner, and with convincing 
accuracy. 
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3 Backlund Transformations 



In ref . , it was established that the double scaled unitary matrix models of refs. |2H| had 
an interpretation as continuum string theories with open string sectors. This was done by 
establishing a direct connection between the solutions of the string equations in equation (H)) 
and solutions of the string equations for those unitary matrix models, which are most naturally 
written in terms of the quantities 5*^, where: 

Si, = Ir',[v^ + v'] - vRk[v^ + v'] , (15) 



as: 



+ hSkHz)] + zv{z) = vC . (16) 



The models are organised by the mKdV hierarchy: 

dv 1 

We define uiz^ and v{z\ such that: 



# = \S'M . (17) 



X^\u,v\ = ^7^'[n] ^ z/r - v{z)n[u] = , 



which implies a specific form for v{z) given a u{z): 



]n'[u] T 

V = - 

Noting the identity [T]: 



v= T • (19) 

/■CM 



= ± z/rx± - 7^[M]x; = {v^ + v')n^[u] - hi[u]n"[u] + ^{n'[u]f - u^v^ , (20) 

we see that if the inverse transformation u = v"^ + v' (known in the integrable literature as 
the Miura map) exists, then this is just our original string equation ([T]); and on substitution 
into equation p9|) gives equation (fTBj) with C = 1/2 ± F. Since the unitary matrix models of 
refs. [571 12H1 were originally derived with C = 0, those models turn out to be identified with the 
case r = 1/2. 

Since the transformation is true for both values of C we find that a given u{z) for a value of F 
gives two functions vc{z) and Vi_c{z), which must be related by: 

vl + v'c = vl_c + v[_c , (21) 
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and hence: 

2C — 1 

— '''' 

On the other hand, the string equation (fTBj) has the symmetry v^c = ~vcj aiid so we obtain a 
result which relates v{z) at one value of C to another for C ± 1: 

2C ± 1 

This implies a similar transformation for functions u{z), relating F and Til. We can make 
this explicit by combining equations (|T9|l . (j20|) . (j23|l and m = f ^ + f ' to give the remarkable 
result: 

3 nz'? - 27^7^" ^ suV w + Au^t^ 

where 7?. = 1Z{uy). 

In the theory of non-linear differential equations (such as those encountered in the integrable 
model context) Backlund transformations are extremely useful, since they convert a given 
solution of the differential equation into a new solution^. Strictly speaking, it looks rather 
like this is not what we have here since our string equation is displayed with F explicitly 
appearing, so Mr and uy±i are solutions of different equations. However, the once-differentiated 
string equation (which is in fact the one which appears naturally in many derivations; see the 
beginning of the previous section, and see below), which is 

^n'" - 2un' -u'n = o , (25) 

does not have an explicit appearance of T. The solutions Mr do know of F, of course, since it 
appears (for example) in their asymptotic expansion 0. 

So we have a genuine Backlund transformation for our system. What is particularly interesting 
is that it is not just any such transformation but a special case of the well-known Backlund 
transformations known to change the number of solitons in the KdV hierarchy. We can make 
this explicit as follows: 

We start with the KdV equations 0, which we rewrite as: 

f ^fiU.Pl-(id»-«-i|)«.Pl. (26) 



^These are actually "auto-Backlund" transformations. 
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where now d=4-- 

ax 

We can search for solutions of the form^ £t(x, tk) = tf^uixtl^). We then find (3k = -I /{2k + 1) 
and that (under appropriate rescahng of tk) equation becomes: 

~u~]^zu' = (^d? ~ud-\?jRk[u] , (27) 

where a prime denotes differentiation with respect to z = xt^*'. Re-arranging we find the 
once-differentiated string equation for a particular k. 

This is another way of making the observation (already noted in ref. jH]) that the string equation 
follows a restriction of the KdV hierarchy to scale invariant solutions. This also fits with 
the recoveryPI of the string equation from the operator relation [P, Q] = Q, where P is the 
generator of scale transformations^. (See the previous section.) 

Now note that the KdV hierarchy admits the following well known Backlund transformation 
relating a solution ui to a solution U2'- 

Wx + yx = ^{w -vY + 2~X , (28) 

where = Ui and i/x = U2- 

This can be derived from the (generalised) Miura map u = v"^ + v' + X (where A is a constant), 
which relates the KdV and mKdV hierarchies, by noting that the mKdV flow hierarchy (see 
equation (fTTj)) is invariant under v — > —v. So we can have two solutions of KdV, Ui and U2, 
arising from the same solution, v, of the mKdV hierarchy. 

We have Ui = = v'^ + v' + X and U2 = Ux = f ^ — f ' + A, which gives upon addition and 
subtraction: 

Wx + yx = 2v'^ + 2A , Wx-yx = 2v' . (29) 

The latter can then be integrated once and substituted into the former to give equation (j2H)). 
This can then be specialized to solutions of the string equation by writing Ui = t^^'^Ui{xt^*') = 
t^^*=//(xt'^*), which gives: 

w = t^'=/i(xt^'=) , y = t'^'^ f2{xt^') , (30) 

®More generally we search for solutions of the form u{x, tk) — t'^''u{xt^''); but we find that consistency forces 
us to choose = 2/3fc. 

^The [P, Q] = Q formalism of ref. l6|_ allows any arbitrary combination of the tkS to be switched on at once. 
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where a prime denotes differentiation witfi respect to z = xt^'=. We liave tlierefore 

/{ + /^ = ^(/i-/2)^ + 2At-^^'=. 



(31) 



So we see that for consistency we must set A = 0, which has an interesting interpretation to be 
discussed below. Hence we have 



So far, it is not clear that this transformation changes T. In order to establish the connection 
between this transformation and the one displayed in equation ()24|) , in which F appears explic- 
itly, we first work perturbatively. Starting with the asymptotic expansions ©, for (say) ui, 
it is easy to show that using equation (jH^ the asymptotic expansion obtained for U2 is of the 
same form, except that T has been replaced by F ± 1. We then expect that they are equivalent 
non-perturbatively, and have checked that this is the case by working numerically on some 
explicit solutions. 

So we have established that the Backlund transformation takes us between solutions with 
asymptotics given in equation © which are for Fs differing from each other by unity. It is 
known |TH] that the Backlund transformations (|28p of KdV increase or decrease the soliton 
number of a given solution by unity. In fact, for a given solution u, the soliton corresponds to 
a bound state of H with (negative) eigenvalue A, and the speed of the soliton is given by —A 
(up to a positive numerical constant). The soliton number of a particular solution can be read 
off by the number of such bound states and their individual velocities are set by the discrete 
spectrum. 

In our case, we have established (so far) that there are no bound states in our solutions, and 
indeed, our restriction of the Backlund transformation forces A = 0. So we deduce from this 
that F, which counts soliton number, must count the number of zero-velocity solitons. This is 
another sign that F is naturally integer, and positive^. We will shortly find further evidence 
for this. 

Now that we have a method for generating new solutions u{z) starting from old ones, we can 
return to our numerical study armed with more powerful tools. 

®It is amusing to note that one can generate non-trivial solutions of the string equation by starting with 
the trivial solution u = 0. This does not give solutions with the "physical" asymptotics given in equation 
which is not a contradiction, since we did not start with such a solution. What one obtains are pole solutions 
with the half-integer values of F similar to those discussed in refs.fTllTS]. 




(32) 
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4 The Case of T = -1 



We established in the last section that given a solution ur, we can generate a solution tir±i- 
In the section before, we reported that it was difficult to solve the string equation numerically 
as r became more negative. Numerical precision was lost rapidly as F went below the value 
—0.9, while positive F is very much under control. We can use our transformation from the 
previous section to surmount this obstacle, by simply solving the equation numerically for u{z) 
for F = e, for e small and positive, and then use our transformation ()24|) to build a solution for 
u{z) with F = — 1 + e. We can simply follow e to vanishingly small values and therefore learn 
the properties of u{z) on the approach to F = — 1. 

We carried this out with very interesting results. See figure El for examples. We were able to 
use this method to generate u{z) for F = — 1 + e where e could easily be taken as small as 
~ 10^^''. The potential well is observed to get more deep and narrow increasingly rapidly as 
e 0. In fact, we observe numerically that the well runs to infinity along the line u{z) = —z, 
becoming infinitely deep and narrow in the limit F = — 1. Most interestingly, note that since 




Figure 3: Numerical solutions to equation for u{z) in the case k = 1. The set of curves on the left 
shows u{z) for a range of negative values of the parameter F, closely approaching the value —1 from 
above. The curve with the lowest well (very close to F = —1) is beginning to show signs of numerical 
error to the right. The curve on the right is the case of F = +1, used for comparison in the text. 

the deepening of the well is accompanied by a narrowing, the system has a chance of preventing 
the appearance of a dangerous bound state. We used the families of u{z) functions found here 
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as potentials for Ti and were indeed able to verify with some confidence (to within acceptable 
error tolerance) that there were no bound states present for any of the potentials. 

So the case F = — 1 is remarkably special, and in fact contains a surprise. At first, it seems that 
it develops a pathology, since as e ^ the well was seen to grow infinitely deep, but at the same 
time it became extremely narrow and, moreover, it moved off to infinity to the right. So the 
solution may in fact be well-behaved. In support of this, it is interesting to note the striking 
fact (again, from studying families of curves such as those displayed in figure Ej) that in the 
limit the F = — 1 curve looks rather like the F = +1 curve! The feature associated to the well 
is confined to a highly localized narrow (and deep) region, which has moved away. Everything 
that remains of the F = — 1 curve at finite z falls (with considerable numerical accuracy) on 
the F = +1 curve. 

This numerical observation suggests the surprising result that the functions Mr=-i and ur=+i 
are actually identical, and that the discrepancy from the piece at infinity actually disappears in 
the limit. We can prove this for all k. The transformation can be used to generate either of 
these functions starting with the case of Mr=o- Putting this function and F = into the equation 
immediately gives the result that Mr=-i = Mr=+i5 since there is no explicit appearance of F in 
the transformation in order to generate the difference between decreasing vs increasing F. 

Note that below F = — 1 there is a transition to completely new behaviour. There, the system 
no longer has a smooth solution with the asymptotics given in equation ©, and there must be 
a singularity at finite z. The origin of the singularity can be seen (for /c = 1) by considering 
the form of the transformation plj) . The denominator of the transformation is 47^^, where, for 
large positive z, we have that 

7^ = u^-2~-^ + ■■■ , (33) 

and so for F < 0, the k = 1 solution must approach the line u = z from below as z ^ +oo. Since 
u{z) started out above this line (at z — > — oo), it must be that it crosses the line somewhere 
in the interior. If this is the case, then somewhere in the interior, the new solution ur-i must 
have a singularity, since the denominator vanishes there. We expect similar arguments to hold 
for the negative F behaviour of other values of k. 

So we have the very natural physical situation singling out F as a positive integer, since any 
solution u{z) with — 1 < F < can be used as a seed for our transformation (jUj) to generate a 
u{z) solution with poles, which have a problematic interpretation. Furthermore, any solution 
u{z) for a positive fractional value of F can also be used, by applying the transformation 
successively, to generate such solutions, and so they are in the same class as the — 1 < F < 
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solutions. The u{z) for positive integer values of F are therefore special, in that they are not 
connected by our transformation to any solutions with poles. 

5 Summary 

So we see that studying the string equation non-perturbatively using a combination of exact and 
numerical methods reveals that the solutions have non-perturbative sensitivity to the values 
of r. In particular, it appears that for F > —1, the spectrum of the model is continuous and 
bounded from below by zero, which is consistent with perturbative expectations. This means, 
for example, that the wavefunction of the extended brane|21] probe[23], the wavefunction ip of 
Ti, thought of as a function ip{Xs) has connected support. 

For F < — 1, the system does not seem to maintain a continuous, bounded spectrum. F = — 1 
is a special value therefore, marking the beginning of a new phase of the model. It is tempting 
to think of this as a phase transition, but since we have every reason to think of F as a discrete 
parameter, this is probably not well motivated. 

We are able to reach all positive integer values of F using our special Backlund transformation, 
and we observe that the cases of F being a negative integer do not really exist, since F = — 1 
is the same as F = +1, so successive application of the transformation will not generate any 
new negative integer cases. So the positive integer F solutions constitute a very special set. 
Meanwhile fractional values of F can be connected by the transformation to solutions with 
poles, and so seem less physical. So while we have not rigourously proven that the positive 
integers are the only physical values selected by the string equation, it is certainly encouraging 
to find that the equation has properties that single out the positive integers as special. 

This connects rather nicely to the rectangular matrix model understanding of the role of F 
mentioned in the first section: If M is "row- like" (F > 0), then the combination MM^ will 
have iV + F eigenvalues, and precisely F of them are zero. These F extra zeros correspond to the 
zero-velocity KdV solitons (see below). On the other hand, if the matrix M is a "column-like" 
(F < 0) rectangular matrix {i.e., more columns than rows) then the combination MM'^ will 
have A^ — |F| eigenvalues. So negative F has no extra zero eigenvalues, and does not correspond 
to having any background branes. In fact it seems to reduce the number of eigenvalues of the 
matrix, but at large N (the limit in which we have defined the continuun model) this makes 
no difference. It would be interesting to determine what the negative F regime may mean 
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physically^. This regime presumably has a plethora of bound states, which would need to given 
a physical interpretation. Perhaps it is simply not possible to do so, and so the development 
of singularities in u{z) for fractional F would be indicative of a non-perturbative sickness, 
reminiscent of that of the even k bosonic models PUI ITTl IT^ . 

The observation that there are F extra zero eigenvalues for positive integer F fits very well 
with the fact that our specialization of the Backlund transformations creates and destroys KdV 
solitons with spectral parameter A = 0. It is clear that the localized D-branes of the minimal 
string theory should be identified with these solitons. This direct connection of D-branes with 
the solitons of the underlying integrable system is intriguing and may well lead to interesting 
new physics in both string theory and in the theory of integrable systems. It deserves further 
investigation. 
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